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If the Higgs mass m# is as low as suggested by present experimental information, the 
Standard Model ground state might not be absolutely stable. We present a detailed 
analysis of the lower bounds on mjj imposed by the requirement that the electroweak 
vacuum be sufficiently long-lived. We perform a complete one-loop calculation of the 
tunnelling probability at zero temperature, and we improve it by means of two-loop 
^ ■ renormalization-group equations. We find that, for ran = 115 GeV, the Higgs potential 
develops an instability below the Planck scale for m t > (166±2) GeV, but the electroweak 
vacuum is sufficiently long-lived for m t < (175 ± 2) GeV. 
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1 Introduction 



If the Higgs boson if sufficiently lighter than the top quark, radiative corrections induced 
by top loops destabilize the electroweak minimum and the Higgs potential of the Standard 
Model (SM) becomes unbounded from below at large field values. The requirement that 
such an unpleasant scenario be avoided, at least up to some scale A characteristic of some 
kind of new physics |TJ, [|, leads to a lower bound on the Higgs mass m# that depends 
on the value of the top quark mass m t , and on A itself. The most recent analyses of this 
bound Q, performed after the discovery of the top quark, led to the conclusion that if the 
Higgs boson was clearly observed at LEP2 (or if m H < 100 GeV), then new physics would 
have to show up well below the Planck scale, Api = 1.2 x 10 19 GeV, in order to restore 
the stability of the electroweak minimum. We now know that win must be larger than 
about 113 GeV [§]]. If m H lays just above this bound, as hinted by direct searches || and 
consistently with electroweak fits ||, absolute stability up to the Planck scale is possible, 
provided m t is close to the lower end of its experimental range [[?[]. For m t around its 
central value, the SM vacuum may not be absolutely stable, but still sufficiently long- 
lived with respect to the age of the universe. Motivated by these observations, we decided 
to reanalyse in detail the lower limits on m# imposed by the condition of (meta-) stability 
of the electroweak minimum. 

We assume that no modifications to the Standard Model occur at scales A smaller than 
the Planck scale. In general, field-theoretical modifications invoked to stabilize the SM 
potential at scales A <C Api, such as supersymmetry, or the introduction of extra scalar 
degrees of freedom, induce computable corrections of order A 2 to the squared Higgs mass, 
thereby forcing A to be of the order of the electroweak scale by naturalness arguments. 
On the other hand, it cannot be a priori excluded that the uncomputable gravitational 
corrections of order A P1 vanish. 

Three different classes of bounds have been discussed in the literature f2[: 

i) absolute stability; 

ii) stability under thermal fluctuations in the hot universe; 

iii) stability under quantum fluctuations at zero temperature. 

The condition of absolute stability is the most stringent one. However, although appealing 
from an aesthetic point of view, this constraint is not demanded by any experimental 
observation: it is conceivable that we live in an unstable vacuum, provided only that it 
is not "too unstable". The condition (ii) — less stringent than (i) and more stringent 
than (iii) — relies on the assumption that the early universe passed through a phase 
of extremely high temperatures (the most stringent bounds are obtained for T ~ Api). 
Although plausible, this is just an assumption; so far, it has been indirectly tested only 
for temperatures up to few MeV. A naive extrapolation of big-bang cosmology by ~ 20 
orders of magnitude in temperature would not only give a bound on the SM Higgs mass; it 
would also exclude various popular unified, supersymmetric or extra-dimensional models, 
because of over-abundance of monopoles, gravitinos, gravitons, respectively. 
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Finally, the requirement of sufficient stability under quantum fluctuations at zero 
temperature gives the less stringent bounds, but does not rest on any cosmological as- 
sumptions. The only cosmological input required is an approximate knowledge of the age 
of the universe Tjy; the bound is formulated by requiring that the probability of quantum 
tunnelling out of the electroweak minimum be sufficiently small when integrated over this 
time interval. In this work we will mainly concentrate on this scenario. 

The probability that the electroweak vacuum has survived quantum fluctuations until 
today is given, in semi-classical approximation, by M 

p « {Tu/RYe~ s ° , (1.1) 

where So is the Euclidean action of the bounce, the solution of the classical field equations 
that interpolates between the false vacuum and the opposite side of the barrier, and R is a 
dimensional factor associated with the characteristic size of the bounce. The main purpose 
of this paper is to reduce the theoretical uncertainties in the above result, performing a 
complete one-loop calculation of the action functional around the bounce configuration 
[ |TT| . As we shall show, this calculation allows us to unambiguously fix the pre-exponential 
factor and the finite corrections at the one-loop level, and also to consistently resum (by 
means of renormalization group equations) the sizeable logarithmic corrections appearing 
in the exponential factor. 

The paper is organized as follows: in Section 2 we shall briefly recall the semi-classical 
result for the tunnelling rate, applied to the case of the SM Higgs. In Section 3 we 
shall discuss the general properties of the one-loop formula, emphasizing the differences 
with the semi-classical one. Section 4 contains all the technical details of the calculation, 
whereas the numerical bounds on the Higgs mass are presented in Section 5. Finally we 
summarize our results in Section 6. 



2 Tree-level computation of the tunnelling rate 

The Standard Model contains a complex scalar doublet <j) with hypercharge —1, 



and tree-level potential 

v{4>) 



(H + iG)/V2 
G~ 



m 2|0|2 + A | |4 = ^_ m 2 R 2 + l _ XH A 



(2.1) 



(2.2) 



where the dots stand for terms that vanish when G, G~ are set to zero. The neutral 
component H is assumed to acquire a non- vanishing expectation value (H) = v. With this 
normalization, v = (G^y^) -1 ^ 2 = 246.2 GeV, and the mass of the single physical degree 
of freedom H is m 2 H = V"(H)\h= v = 2Xv 2 . As is well known, for H ^> v the quantum 
corrections to V(H) can be reabsorbed in the running coupling A(/i), renormalized at a 
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scale pL ~ H. To good accuracy, V(H » v) = \\{H)H 4 and the instability occurs if, for 
some value of H, X(H) becomes negative. Since, for mn larger than 100 GeV, this occurs 
at scales larger than 10 5 GeV, we shall neglect the quadratic term m 2 H 2 /2 throughout 
the paper. 

In general, the bounce || is a solution H = h(r) of the Euclidean equations of motion 
that depends only on the radial coordinate r 2 = x^x^. 

- dltdll h + nh) = ~-~ + V{h)=0, (2.3) 
and satisfies the boundary conditions 

h'(0) = , h(oo) = v -> . (2.4) 
We can perform a tree-level computation of the tunnelling rate with A < 0. This leads to 

, / s / 2 2R _ r , , 87r 2 ,„ . 

h(r) = iw^TW • s " lh] = m ' (2 ' 5) 

where R is an arbitrary scale. At first sight, the approximation of taking V(h) = Xh 4 /4: 
may appear rather odd, since the unstable vacuum configuration h = corresponds to 
the maximum of the potential. However, this is not a problem within quantum field 
theory, since the tunnelling configuration requires a non-zero kinetic energy (the bounce 
is not a constant field configuration) and is therefore suppressed even in the absence of a 
potential barrier ||. The SM potential is eventually stabilized by unknown new physics 
around Api: because of this uncertainty, we cannot really predict what will happen after 
tunnelling has taken place. Nevertheless, a computation of the tunnelling rate can still 
be performed ||. 

The arbitrary parameter R appears in the expression of the bounce since, because of 
our approximations, the potential is scale-invariant: at this level, there is an infinite set 
of bounces of different sizes that lead to the same action. 

Substituting the bounce action Q2.5Q in Eq. ( |1.1|) , the condition p < 1 for a universe 
about 10 10 years old is equivalent to A > —0.065/(1 — 0.011ni?u), i.e. A cannot take 
too large negative values. The bound on A can be translated into a lower bound on 
m H taking into account the renormalization-group evolution (RGE) of A(/x) (see Fig. |l|). 
At this stage, however, there is clearly a large theoretical ambiguity due to the scale 
dependence. Which values of R and of the RGE scale \i should one use? As we shall 
see in the next section, both ambiguities are solved by performing a complete one-loop 
calculation of the tunnelling rate. 

3 Quantum corrections to the tunnelling rate 

The procedure to compute one-loop corrections to tunnelling rates in quantum field theory 
has been described by Callan and Coleman |I0| , following the work of Langer []ITJ in 
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statistical physics. At this level of accuracy, the tunnelling probability per unit four- 
dimensional volume V can be written as 



P_ 
V 



-Si[h] 



S 2 [h] 
4vr 2 



SDet' S ' [ft] 



SDetS '[0] 



-1/2 



Sol 



(3.1) 



where h still denotes the tree- level bounce and So (Si) is the tree-level (one-loop) action 
functional. Here h = indicates the false (electroweak) vacuum; we assume that the 
potential has been shifted so that Si[0] = 0; S ' denotes double functional differentiation 
of So with respect to the various fields; Det is the functional determinant, and SDet = Det 
or SDet = 1/ Det 2 , depending on whether it acts on boson or fermion fields. 

When evaluated with a constant field configuration, Si is simply given by the usual one- 
loop effective potential. Computing Si [h] is a much harder task because the bounce (|2.5| ) 
is not a constant field configuration. Furthermore, unlike the constant-field case, there 
are quantum fluctuations that correspond to translations of the bubble. The 'prime' on 
SDet' S ' [/i] in Eq. ( |3.1| ) indicates that these fluctuations, corresponding to zero modes, 
have been explicitly removed from the functional determinant. In this way the result 
acquires a dimensional factor that will be compensated by the integration over the volume 
of the universe. 

Fortunately, in order to compute the one-loop corrections to the tunnelling rate we 
do not need to find the field configuration hi that extremizes the full one-loop action: 
we only need to compute Si[/i], where h is the field configuration that extremizes So and 
has the simple form in Eq. ( |2.5| ). The difference between S\[h] and Si[/ii] is a two- loop 
correction. 

In our case, the main effect of quantum fluctuations is the breaking of scale invariance 
of the tree-level potential. As we shall show explicitly, this implies that bounces with 
different R, which have the same action at the semi-classical level, turn out to have 
a one-loop action roughly given by Si[h] ~ 87r 2 /(3|A(l/i?)|). At the same time, for 
these configurations the dimensional factor due to the zero eigenvalues turns out to be of 
0(R~ A ). In this way the two scale ambiguities of the semi-classical result are completely 
resolved. Indeed, the complete result for the tunnelling probability at one loop can be 
written as 

8tt 2 



V 



Vu 

max — ■ exp 
R i? 4 F 



'3 | AMI 



-AS 



(3.2) 



where 



AS 



-3(5 + 6L) 

_2gl + gl L 
32A 2 1 



9t 
6IA 



(13 + 12L) 



2L) + f h (\) 



9t 
6A 2 



(5 + 6L) 



2g 2 2 



9 2 z 



121 A| 



(7 + 6L) 



(3.3) 



Here Vj 



ln(i?/ie 7E /2), g t is the top Yukawa coupling, and g%, g% 
the weak gauge coupling constants, defined at tree level by m t = g t v/y/2, r 
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V gy are 
S 2 V/4, 

g 2 z v 2 /4. All couplings are renormalized at the RGE scale /i. The numerical values 
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Figure 1: Running of the quartic Higgs coupling for mn = 115 GeV and m t = 
165, 170, 175, 180 and 185 GeV [a s (m z ) = 0.118]. Absolute stability [A > 0] is still possible 
if m t < 166 GeV. The hatched region indicates the metastability bound. 



of the functions f g , f t are plotted in Fig. ||, whereas fh is given in Eq. (|4.61|) . The L terms 



cancel the fi dependence of A in the leading semi-classical term. If one chooses a value of 
/i ~ 1/-R, such that L ~ 0, the typical correction to the action is of 0(gf/X 2 ) ~ 10, to be 
compared with the leading term of order 100. 

In previous analyses (see e.g. Ref. a full one-loop computation of the tunnelling 
rate was never performed, and the semi-classical result was improved by considering only 
quantum corrections to the effective potential, or to the running of A. This procedure 
leads to a correct estimate of the leading logarithmic corrections to the action, but the 
finite terms of the calculation are not under control. Within this approximation the use 
of two-loop RGE equations does not improve the accuracy of the calculation. On the 
other hand, a consistent implementation of two- loop RGE equations for A(/i) is possible 
starting from Eq. (|3.2|). 

In Fig. |l] we plot the evolution of A(/i) as obtained by integrating the two- loop RGE 
equations of A, the top Yukawa coupling g t and the three gauge couplings g, L for 
rriH = 115 GeV and some reference values of the pole top mass m t .f\ For comparison we 
also show the lower bound on A derived from Eq. ( |3.2| ), imposing the condition p < 1 and 
assuming Vjj = (10 10 yr) 4 . As can be noticed, the evolution of A crosses the metastability 

lr The initial values of A and and g t have been related to the values of run and mt using the matching 
conditions given in [Q and |15| , respectively. The discussion about the uncertainties involved in this 
estimate of X(fi) is posponed to Section 5. 
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Figure 2: Numerical results for the subtracted part of the correction to the action from top 
loops [ft(gt/\M)> S0 Hd curve] and gauge boson loops [fg(g 2 /\M), dashed curve}. Relevant 
values of g 2 /\X\ are 4 8. 



bound (i.e. the tunnelling rate becomes too high) for values of /i much larger than the 
electroweak scale. This implies that our approximation of neglecting the 0(v 2 ) quadratic 
term in the tree-level potential is very good, since the critical bounces are those with a size 
much smaller than 1/v. It is also important to notice how the lower bound on A increases 
as a function of the RGE scale (or oil/ R). This effect is due to the pre-exponential factor 
in Eq. ( |3.2| ), scaling like R~ A , which we have been able to determine from the one-loop 
computation. It is important to notice that, for the experimentally interesting values of 
nin and m t , the tunnelling rate is dominated by bubbles with 1/R about two orders of 
magnitude below Api, as can be seen in Fig. [I] or, more clearly, in Fig. |3]. Therefore the 
metastability bound on tuh does not depend on the unknown physics around Api. 

4 Explicit computation of the one-loop action 

4.1 General strategy 

The central point in the computation of the tunnelling probability, Eq. ( |3. l| ) , is the eval- 
uation of ratios of functional determinants: 



in the various sectors of the theory. This requires solving eigenvalue equations of the type 



DetS^h] 
DetS#[0] 



(4.1) 



SZ[h}^ = i-d 2 



+ W(r)]iJ> = 



(4.2) 
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Figure 3: Contribution to the tunnelling rate (in arbitrary units) from bubbles of different 
R, as a function of 1/R, for mjj = 115 GeV and m t = 175 GeV. 



where ifj generically denotes scalar, fermion or gauge fieds. In order to perform such a 
calculation, one should i) choose a suitable eigenfunction basis; ii) define the renormaliza- 



tion procedure. In both respects, our approach will be similar to that of Ref. |16[ , where 
sphaleron computations within the SM have been performed. Details will be given in the 
next subsections; here we just sketch the main points of our strategy. 

Because of the four-dimensional spherical symmetry of the bounce, the 'interaction 
term' W in ( [1.2[ ) depends only on the radial coordinate r = (x^x^) 1 ^ 2 . For this reason 
it is convenient to decompose the various fields in eigenstates of the four-dimensional 
angular momentum operator L 2 , and to write the Laplace operator as 

d 2 ^ r\ J 2 

9 2 = 9A = ^ + -^-^-V L . (4.3) 
ar z r dr r z 

In the case of scalar fields, the eigenfunctions of L 2 are the four-dimensional spherical 
harmonics Yj(6) (where 9 collectively denotes the 3 polar angles) with eigenvalues 4j(j+l), 
and degeneracy (2j + 1) 2 , where j takes integer and semi- integer values (see Appendix |AT). 
After this decomposition, we have 

logDet^' = $>gDet(Sa- , (4-4) 

j 

where (S'^j is the restriction of Sq to the subspace spanned by eigenfunctions with angular 
momentum j. As we shall see, the situation is slightly more complicated for fermion and 
vector fields. Their expansion in spinor and vectorial hyperspherical harmonic functions 
will be obtained starting from the Yj(8). 

A further simplification arises from the fact that, in order to compute the ratio in 
Eq. ( fO]) , it is only necessary to solve Eq. ([OP for A = 0. Indeed, we can use the 
result lITi 



Det [-V L + W{r)\ detti 



Pi = " r „ v /J = lim 7 V ' , (4.5) 
Hj Det[-V L ] ™ detu J (r) V ; 
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where u ] w (r) [u (r)] are eigenf unctions, regular at r = 0, of — Vl + W(r) [Vl] with zero 
eigenvalues: 



[-Vl + W(r)\ u J w (r) = 0, V L ui(r) = . (4.6) 
The symbol 'det' in Eq. fl4.5|) stands for the ordinary determinant over residual (spinorial, 



16 and the next subsections for 



gauge group, etc.) indices of these solutions (see e.g. 
more details). 

The one-loop action is affected by the usual ultraviolet divergences of renormalizable 
quantum field theories. As a consequence, the sum over j in Eq. (|4.4j ) is not convergent 
(the ultraviolet behaviour being encoded in the behaviour of the determinants for j — > oo), 
and the usual renormalization procedure is needed. The expression 



Si = S + i In SDet S^[h] - i In SDet S£[0] 



So + AS 



(4.7) 



can be made finite by adding an appropriate set of local counterterms, which lead to a 
redefinition of the bare couplings in So- For example, in the MS scheme (dimensional 
regularization with MS subtraction) the renormalized one-loop action can be written as 



S\ — S { 



MS 



AS- (AS) 



pole 



(4.8) 



where S^ is the lowest-order action expressed in terms of the renormalized MS couplings, 
and (AS') pole is the divergent part of AS" defined according to the MS renormalization 
prescription. Computing the full determinant in d = 4 — 2e dimensions would be an 
extremely difficult task. However, this is not necessary. In fact, the divergent terms 
are all contained in AS 1 ^, defined as the expansion of AS up to second order in the 
'interaction' W: 



AS [2] = - [in SDet (-d 2 + w) - In SDet (-9 



STr 





- -STr 




4 



0(W 2 ) 

d^wi-dY'w 



(4.9) 



where STr = Tr or STr = — 2Tr depending on whether it acts on boson or fermion fields. 
In other words, the difference AS" — AS^ is ultraviolet-finite. Equation (|4.8|) can be 
rewritten as 



Si 



qMS 
°0 



AS - AS [2] 



AS [2] - (AS) 



pole 



(4.10) 



where the two terms in square brackets are separately finite. The advantage of this last 
expression is that AS^ can be computed either as a (divergent) sum of terms corre- 
sponding to different values of the angular momentum, which gives a finite result when 
subtracted from AS, or by standard diagrammatic techniques in 4 — 2e dimensions. In 
the next subsections we shall show how this procedure is implemented in practice. 

The final result is expressed in terms of the renormalized parameters A, gt and gi, 
whose definitions depend on the renormalization scheme. Of course, the scheme depen- 
dence disappears once these couplings are re-expressed in terms of physical observables, 
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such as Higgs and top pole masses. In practice, however, in the case of the gauge cou- 
plings it turns out to be more convenient to directly use the MS definitions, since these 
parameters are accurately determined by fitting multiple observables. 



4.2 Higgs fluctuations 

The relative corrections to the action due to fluctuations of the Higgs field are generally 
small because the Higgs coupling A is small; it is however important to consider them, 
because they include the special zero-modes (eigenfunctions with zero eigenvalue) corre- 
sponding to translations (j = 1/2) and field dilatations (j = 0) of the bounce, as well 
as the unique negative eigenvalue corresponding to space dilatations of the bounce. It is 
precisely the existence of this negative eigenvalue that makes the false vacuum unstable. 
In this case the interaction term is simply given by 

94 R 2 

W(r) = V"(h) = • (4.H) 

It does not depend on any coupling constant, as expected, since the leading contribution 
to the action is proportional to 1/A. 

The solutions of the 'free equation' Vl^oM = 0; regular in r = 0, are immediately 
found to be proportional to r 2 - 7 (see Eq. ( |4.3| )). In order to exploit the result in Eq. ( f4.5| ), 
we must compute the ratio 

pAt) = 4tt • ^ 

u J (r) 

The functions Pj{r) obey the differential equations 

pj(r) + ^f^(r) = W{r) Pl (r) , (4.13) 
which can be solved analytically: 

P^r) = {r l + R 2 ? > ( 4 - 14 ) 



where 



2j-l , j(2j - 1) r/M _, 

a; = ; bj = — r . 4.15 

3 J + l 3 (j + l)(2j + 3) 1 J 



According to Eq. (|4.5|) , we have therefore 



Pet [-V L + W(r)} 

Pi = ? ; = hm pAr) = b« . 4.16 

Hj Det [-V L ] ) J v > 

Let us neglect for the moment the contributions of j = and j = 1/2, for which pj = 0. 
Taking into account the multiplicities of the sub-determinants, we have 



i>l/2 j>l/2 I 



2j Af 



(4.17) 
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As discussed in the previous section, we regularize this expression by subtracting from it 
the terms obtained by solving Eq. ( [4.1 3|) perturbatively in W. Replacing W with eW, 
we define the coefficients hij as Pj(e) = 1 + eh±j + e 2 ti2j + 0(e 3 ). These coefficients can 
be easily determined numerically. The subtracted series is rapidly converging, and in 
practice the inclusion of the first ten terms already provides an excellent approximation 
of the full result. We find 



\ J2 (2j + l) 2 lnp,-i£(2j 

Z j>l/2 Z j>0 



h 



^i] = 12.6 



(4.18) 



We now turn to a discussion of zero eigenvalues. The four zero-modes in the j = 1/2 
sector correspond to translations of the bounce. They can be converted into a volume 
factor following the procedure illustrated in [1C], which amounts to replacing pu 2 with 



, Det/(gff[ft]) 1/2 Pet [e+QSff[fe]) 1/2 ]/e 

Pl ' 2 Det(^[0]) 1/2 ™ Det(5ff[0]) 1/2 



(4.19) 



and multiplying the expression for the tunnelling probability per unit volume by a factor 



So[h]/2ir for each of the four translation zero modes. This is the origin of the factor 
S 2 [h]/(2n) 2 in Eq. (|3.1|), which, in our final result Eq. (|3.3|), is included in f h . The 
elimination of the four vanishing eigenvalues from pin provides the dimensional factor 
R~ A in Eq. ( |3.2| ). In fact, after numerical integration of the corresponding differential 
equation, we find 

p' 1/2 = 0M1R 2 . (4.20) 

As can be see from Eq. fl4.15| ), there is another zero eigenvalue in the j = sector. It 
arises from the scale invariance of the tree-level potential: bubbles with different field value 
~ 1/R have the same tree-level action. Scale invariance is broken by quantum corrections, 
which shift the zero eigenvalue by an amount of 0[g 2 /(AttR) 2 }. We can therefore compute 
Po by removing the zero eigenvalue from the determinant as in Eq. ( [4.19|) , and replacing 
the zero eigenvalue with a quantity of 0[gf /(AirR) 2 ] (this small correction can in principle 
be extracted from the R dependence of the one-loop bounce action [|18|). This gives 
-R 2 , and 



ln |Po| 



In 



Po- 



91 



-2.5 



(4.21) 



' (47T.R) 2 

for p' Q = —R 2 and gt — 1. Note that p' Q is negative, as it should be because of the 
instability of the electroweak vacuum |l(J . 

Combining the results in Eqs. ( [4.18| )-( f4.21| ) we finally obtain 



= \(AS'-AlnR) - AS [2] 



Higgs 



4± 1 



(4.22) 



where the quoted error is entirely ascribable to the uncertainty in Eq. ( [4.211 ). In Sec- 
tion O, will be combined with analogous contributions generated by Goldstone 



boson fluctuations to yield the constant fh in Eq. 
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Renormalization 



Following the general strategy discussed in Section 3~T , we now proceed to evaluate AS 1 ' 2 ' 
in the MS scheme. The explicit <i-dimensional expression of AS®, defined in Eq. Q4.9|) , is 

3A/i 4 " d . r .,, , 9AV {4 " d) 



[AS® 
where 



Higgs 



Ai J d d x h\x) 



J d d xd d y h 2 (x)A 2 (x-y)h 2 (y) , 

(4.23) 



Ai 



A 2 (x-y) = J 



d d k 1 
(2 7 r) d k 2 ' 

d d q 
W) 



iq(x-y) 



d d k 



(2%) d k 2 (k + q) 2 ' 



(4.24) 
(4.25) 



and n is the usual mass scale that is introduced in dimensional regularization in order to 
keep the action dimensionless. In dimensional regularization, Ai = 0, and there is thus 
no contribution to AS® from terms of order 0(W). This is also true for the 0(W) terms 
in the fermion and gauge sectors, as a consequence of our approximation of neglecting 
mass terms in the tree-level potential. The second term in Eq. ( |4.23| ) can be written as 



AS® 



9A 2 //" 



Higgs 



d d q 
(2tt) c 



h 2 (q 2 ) B (q 2 



where h 2 (q 2 ) is the <i-dimensional Fourier transform of h 2 (x) and 



A-d 



d d k 



(27r) rf k 2 (k + q) 2 (4tt) 



1 , IL A 

- + 2 + In ^ 

e q z 



Oil) 



(4.26) 



(4.27) 



Here 1/e = 2j{d — 4) — 7^ + ln(47r) denotes the divergent part, to be subtracted according 
to the MS prescription. We therefore have 

1 9AV~ d 
i 64vr 2 



(A5')f iggs 



d d q 



(2vr) d 



/i 2 (g 2 



(4.28) 



Since 



d d q 
(2tt) c 



/i 2 (g 2 



d d z ^ 4 (x) 



(4.29) 



the divergent part in Q4.26 ) is immediately recognized to have the same structure as the 
bare action. 

After subtraction of the 1/e pole, the integral in Eq. ( |4.26| ) can be safely performed 
in four dimensions (see Appendix jB|), leading to 



AS [2] - (AS) 



pole 



9A S 



d 4 q 

Hi ggs 64vr 2 J (2vr)^ 



h 2 {q 2 



In 



q 2 



-3L 



5 
2 



(4.30) 



where L = ln(i?/xe 7E /2). Combining the above result with Eq. Q4.22 ), we obtain the full 
one-loop correction to the action due to fluctuations of the Higgs field, renormalized in 
the MS scheme. 
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4.3 Top fluctuations 



The fermionic determinant due to fluctuations of the top-quark field assumes a form 
similar to the bosonic one, if expressed in terms of the squared Dirac operator: 



9t ,2 



9t 



_q2 + « ^ + jn 



-d 2 + W 



(4.31) 



(4.32) 



x 2 " ' V2 J 
where g t is the top Yukawa coupling. Indeed, we can write 

A 1 S top = -^[lnDet (flip) -I: 

where N c = 3 is the number of colours. 

The only difference with respect to the Higgs case is the structure of W. Since W does 
not commute with L 2 because of the (j)h term, we cannot decompose the eigenfunctions 
as products of scalar spherical harmonics times constant Dirac spinors. However, the 
spherical symmetry of the bounce, which implies 



<jHh(r) = $h'(r) , x^ = x^/r 



(4.33) 



considerably simplifies the problem. Following Ref. |L9| , we decompose the fermionic field 



as 



(4.34) 



[au(r)^u(0) + a 2 j{r)%l) 2 j 

J=±3/2, ±5/2,... 

where atij(r) are scalar functions and ipij{9) are spinors written in terms of appropriate 
hyperspherical harmonic functions of multiplicity J 2 — 1/4 (see Appendix ^). In this basis 
the equation fllpip = becomes 



<r 2 



dr 2 + r 2 



+ f/i 2 (r) 



_ d 2 , J(J+1) , 9? h 2, 
dr 2 ~+~ r 2 ~+~ 2 fl \ 



Oi2j{r) 







(4.35) 



In the limit gt = the two components of Eq. (|4.35|) are decoupled; for J > their 
solutions, regular in r = 0, are given by a\j oc r J and a^j oc r J+1 . Since Eq. fl4.35| ) is 
invariant under the exchange J <-> — J, we shall consider in the following only the case 
J > 0, and include an extra factor of 2 in the multiplicity of the solutions. Equation Q4.35Q 
can be cast in the form 



a i o / y * u2 i y* i.i 
Pu + 2 -Pu = Y Pu + ^hp 2J r 

p'L + ^p'v = ^ PaJ + 



gt _i 

hpijr 



where pu(r) = aij(r)/a^j(r). The system Eq. 
with initial conditions pfj(0) = Sik- According to Ref. 



AS, 



top 



-3 E 

J=TN+l/2 



J — — ] In i lim det 



has two solutions pfj(r) 
we can finally write 

p\A r ) Pij(r) 
p\j{r) Ph(r) 



(4.36) 



1,2 



(4.37) 
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The sub-determinants at fixed J can be computed by means of numerical methods. As 
usual, the sum in Eq. ( |4.37|) is divergent and we regularize it by subtracting the first two 
powers in W . Since g t appears always multiplied by h(r), the final result depends only 
on the ratio <?f/|A|: 



AS - AS [2] 



top 




(4.38) 



and is plotted in Fig. |] in the range of interest. 

The evaluation of AS^ in dimensional regularization is very similar to the Higgs case. 
We find 



[AS® 



d _ Nc 
top 4 



d d q 
(2tt) c 



{tif-'W)? + 2gh 2 ~h 2 (q 2 )} B (q 2 ) , (4.39) 



which, after subtraction of the divergent terms, leads to 



AS® - (AS) 



pole 



9t 



top 



6 A 



(5 + 6L) 



_9t_ 
6|A| 



(13 + 12L) 



(4.40) 



4.4 Gauge and Goldstone fluctuations 

In order to evaluate quantum fluctuations in the sector of gauge fields we need to mod- 
ify the bare action introducing an appropriate gauge fixing. Both the correction to the 
potential and that to the kinetic term are gauge-dependent, but they combine to give 
gauge-independent physical quantities pOfl . For example, at one loop, the divergent cor- 
rections change the Higgs action into 



S 



d x 



;i + z T )|^0| 2 + (i + z y )A|0r 



(4.41) 



and the gauge dependence of Zt and Zy cancels in the combination Zy — IZt that 
determines the RGE equation for A. By requiring that the bounce action be an extremum 
under the transformation h(r) — > /i(£r), one finds that St = —2Sy, where St (Sy) is the 
'kinetic part' ('potential part') of the tree-level bounce action S = Sr + Sy. Therefore, the 
divergent corrections to the bounce action are gauge- independent. More generally, it was 
shown by Nielsen |2T[] that the value of the effective potential at its minimum is gauge- 
invariant. Similarly, the value of the effective action is also invariant, when evaluated over 
a field configuration that extremizes it. It follows that the full corrections to the bounce 
action are gauge-independent. 

The functional obtained after gauge fixing depends on gauge fields A 1 ^, Goldstone 
bosons G b and Faddeev-Popov ghosts if. We are only interested in the second functional 
derivative of the action evaluated at A ^ = G b = if = and h ^ 0. This considerably 
simplifies the problem, since we can ignore the non-Abelian part of gauge interactions. 
In the following, we will denote with a suffix A the corrections to the action for a generic 
abelian gauge field A^, with coupling g. We can also ignore the fluctuations of the 
electromagnetic field: they vanish in the difference AS'f/i] — AS^O] because the photon 
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field is not coupled to h. The problem is thus equivalent to the case of three independent 
Abelian fields: the two Ws and the Z. 

Introducing a 't Hooft-Feynman gauge-fixing term, which has proved to be the most 
convenient choice for our calculation, the relevant part of the bare action can be written 

as 



d x 



4 ^ 



1 



(OA - ghGf + \D<P\ 2 + V{<f>) + v*[-d 2 + M(G, h)]r) 



(4.42) 



where M(0, /i) = g h /4. This choice of the gauge fixing term has the advantage that 
all terms of the type A^d^G and d^d u A^, generated by \D(f>\ 2 and F 2 , respectively, are 
eliminated from the equations of motion. 

The ghost fields can be treated separately since the the second derivative of the action 
is diagonal with respect to them. On the contrary, we cannot separate A^ and G, whose 
zero-eigenvalue equation is given by 



-d%„ + \g 2 h\r)5 i 
-gx u h'(r) 



-gx^h'ir) 
-d 2 + (W + \) h 2 ( 



A u (x) 
G(x) 



. 



(4.43) 



In order to solve Eq. (|4.43|) we need a suitable decomposition of A^ in terms of spherical 
harmonics. A convenient choice is provided by 



E 

J=0,i,l, 



aij(r)x^ + 



a2 ^ rd + 
2[J'(J + 1)] 1/2 " 



i=3,4 



Yj{0), (4.44) 



where the a^(r) are scalar functions and Vj 1 ' are two generic orthogonal vectors. In this 
basis the operator L 2 is diagonal with respect to the index j and, at fixed j, mixes only 
the first two components, a\j and a 2 j. Since 9 M F — 0, the case j = deserves special 
attention. For this reason, we discuss separately the two cases j > and j = 0. 



Sub-determinants with j > 

As shown in Appendix |A|, for j ^ the zero-eigenvalue equation is 



-V 



L 1 j.2 1 4 

-gh'(r) 



Vi(i+i) 
-V L - 4r + \g 2 h 2 (r) 



-gh'{r) 


-v L + a g 2 + \)h 2 (i 



-V L + -g 2 h 2 (r] 



aij{r) 







3,4 



ay(r) " 




a 2j (r) 


= 0, 


Gj{r) _ 






(4.45) 




(4.46) 



where Gj(r) are the components of a standard decomposition of G{x) in spherical har- 
monics. 

Interestingly, the components a 3 j and a^, which are decoupled from the Goldstone 
boson sector, obey exactly the same equation as the two components of the complex 
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ghost field. The latter contribute to the action with an opposite sign; therefore these two 
contributions cancel against each other for j ^ 0. 

Equation (|4.45|) can be further simplified by means of a rotation in the (ay, a 2 j) space, 



a 2j 



1 



VJ -vT 



VT 



VI 



aij 
a 2j 



which diagonalizes it in the gauge sector. In this new basis we have 



-V L+ + 



(4.47) 



-V L _ + \g 2 h 2 



-gtiJj/(2j + l] 



-V L+ + \g 2 h 2 gh'y/V + l)/(2j + 1) 

-gh'y/j/Qj + 1) gh'y/U + l)/(2j + 1) -V L + {\g 2 + \)h 2 











" aij ' 






a 2j 


= 




. Gj _ 





(4.48) 

where Vl ± denote the usual d 2 operator at fixed angular momentum, with eigenvalues 
shited by j — > j ± 1/2. The free equation {h — > 0) is now diagonal and, similarly to the 
top and Higgs cases, can be easily solved. We finally obtain 



" i ' 
Pij + — — Pij 



P2, 



Pi, 



5 + 4j , 

3 + 4j , 
— —Psj 



-g h p^-rp^gti^^j 

1 2r2 1 ./ / J + 1 

ighp^ + ^ghJ — 



\g 2 h 2 + Xh 2 ^j p^-^gh'^j— 



j 



rp 2j gh\ 



I j + 1 
2j + l 



where, as in the previous cases, we have defined 



Pij 



a -' 



1,2; 



P3j 



Gj_ 
Gj 



(4.49) 
(4.50) 



In complete analogy with the top case, the above system has three solutions pjj(r) (k 

} >.i 



1,2,3), with initial conditions pf 5 (0) = 8^, and the resulting correction to the action can 
be written as 



A 1 A 



-£(2j + l) 2 ln^ hmdet 

z j>0 



p\ 3 {r) p y (r) p%{r) 

p\ 3 { r ) P 2 2j( r ) P%(r) 
Ph(r) pUr) pUr) 



(4.51) 



All the determinants in Eq. ( |4.51| ) are different from zero and, as usual, their sum must 
be regularized by the subtraction of AS^ . 



Sub-determinant with j = 

For j = only the x M component of is different from zero [^4|f =0 ' ) = OioM^m] anc ^ the 
eigenvalue equation is 



' dr 2 



3 d , _3_ 

T r 2 

-gh'(r) 



'H + ^ + l9 2 h 2 (r) -gh\r 

_ d?_ 

dr 2 r dr 



3d +a g i + \)h 2 (r) 



aio{r) 
Go(r) 



. (4.52) 
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This leads to 



// 
Pi 



P2 



5 , 

-Pi 

r 

3 , 

-Pi 

r 



4 



g h pi - -p 2 gti 
r 

g 2 h 2 + Xh 2 ) p 2 - r Pl gti 



(4.53) 



where, as usual, p% = aio/a5 and p 2 = Gq/Gq. The absence of transverse components 
in A^ =0 ^ implies that the j = sub-determinant of the ghost fields is not cancelled; it is 
obtained as the large-r limit of the function p v {r), where 



Pn 



3 , 

-Pr, 



2 1,2 

g h p v 



r ' 4 

Employing the usual notation for the two solutions of the system Q4.52j ), we can write 



(4.54) 



A Q A 
^3=0 



- In < lim det 

2 r-»oo 



p\{r) p\(r) 
P2O) PlO) 



In lim Pn(r) . 



(4.55) 



The 2x2 determinant in Eq. ( |4.55| ) contains a vanishing eigenvalue. This is present also in 
the limit g = 0, and originates from the global symmetry of the action, which persists even 
after the introduction of the gauge-fixing term. This corresponds to the possibility, for 
the false vacuum, to tunnel in different directions with equal probabilities. As discussed 
TBI, E 



in 



, the fluctuations corresponding to global rotations can be converted into an 
integral over the group volume. These zero-mode fluctuations probe the non-abelian 
structure of the gauge group, so that we can no longer compute separately the Z, W ± 
contributions. Considering the full gauge group, we obtain 



16tt' 



— / d A x hf 
2tt 



3/2 



-AS' 

gauge 



R' 



f (2) 
J h 



-AS" 

gauge 



(4.56) 



where the factor 167r 2 is the volume of the broken SU(2) group. The resulting integral of h 2 
is logarithmically infrared divergent, as a consequence of our approximation of neglecting 
the mass term in the Higgs potential, which would have acted as an infrared regulator. 
For this reason, we cut off the integral in Eq. ( ^4.56| ) at r = 1/v. The result, for Rv 1, 
is 

"217/3^7/3 ! 

m — 

Rv 



Jh 



(2) 



2 



A 



+ 0{l/lnRv) 



(4.57) 



This vanishing eigenvalue is the only aspect of our calculation that is sensitive to the 
infrared behaviour of the potential, and therefore beyond the control of our approximation. 
However, the sensitivity to the infrared cut-off is mild and does not induce an appreciable 
numerical uncertainty. 

Once the vanishing eigenvalues have been removed from the j = sub-determinant, 
the latter acquires a dimensional factor that is compensated by the factor R 3 in Eq. ( |4.56| ) . 
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Figure 4: Metastability region of the Standard Model vacuum in the (m#,m t ) plane, for 
&s( m z) — 0.118 (solid curves). Dashed and dot-dashed curves are obtained for a s (mz) = 
0.118±0.002. The shaded area indicates the experimental range for m t . Sub-leading effects 
could shift the bounds by ±2 GeV in m t . 



Final result for the gauge sector 

Summing the regularized gauge and Goldstone corrections, we finally obtain 



AS. 



where 



and 



f (3) 

Jh 



gauge 



lim 

9-0 



2/ 9 



|A|, 



+ f 9 



+ fr + 3/, 



h > 



(AS'-]nR)+AS. 



AS [2] 



Mm) = 



\nR) + AS j>0 - AS [2] ] 



(3) 
h ) 



(4.58) 

(4.59) 
(4.60) 



so that f g (0) = 0. The function f g , obtained by means of numerical integration, is plotted 
in Fig. 0. The constants , can finally be combined with f^ in Eq. (|4.22|) and with 
the translation prefactors Sq/{2k) 2 to yield 



, _ f (i) f (2) o f (3) , 16tt 2 

Jh - Jh + Jh + 6 Jh - 111 



lnlAI -8± 1 . 



(4.61) 



This numerical estimate has been obtained with Rv = 10 14 , but it is very mildly sensitive 
to the value of Rv. 
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In analogy with the previous cases, dimensional regularization leads to 



AS^ 



d d q 
{2ix) d 



S-*[h?(q*)]' - 2g¥kV) }B„( q 



which, after subtraction of the divergent terms, reduces to 

A 5 + 6L g 2 7 + 6L g A l + 2L 



AS [2] - (AS) 



pole 



|A| 12 A 2 32 



(4.62) 
(4.63) 



Results 



Using Eq. ( ]3.2|) , the condition p < p max can be translated into an upper bound on |A(/x)|: 



\X{2e^ E /R)\ < 



8tt 2 /3 



HVu/R^-AS-lnp^ 

0.065 



[\n(Rv) + MTft/Tu)] /100 - [AS + \np max ] /400 



(5.1) 



where AS is computed at L — ln(i?/ie 7j5 /2) = 0, Vu — T{j and Ty = 10 yr. Equa- 
tion ( |5.1|) must be satisfied for all values of R. From this inequality, a lower limit on ttlh 
can be obtained by integrating the two- loop RGE equations [13| for A, g t and the three 
gauge couplings The initial values of A and g t at \i — v are related to the values of m# 
and mt by the matching conditions given in [14] and [0, respectively. As in most recent 
analyses of the Higgs potential [H], we include the finite two-loop QCD correction in the 
relation between gt{mt) and the top pole mass m t . The latter is formally a higher-order 
contribution, and can be used to estimate the theoretical uncertainty in the determination 
of A(/i) at high scales. 

Stability and metastability bounds in the (m//,m t ) plane are shown in Fig. [|, which 
has been obtained with p max = 1. The metastability condition can be approximated as 



m H (GeV) > 117 + 2.9 K(GeV) - (175 ±2)] -2.5 

while the absolute stability bound is given by 

m H (GeV) > 133 + 2.0 K(GeV) - (175 ± 2)] - 1.6 



a s (m z ) - 0.118 
0.002 



0.1 In 



10 10 yr 



a s {m z ) -Q.lli 
0.002 



(5.2) 



(5.3) 



The numerical impact of AS is rather weak, provided the renormalization scale is chosen 
such that L = 0: for m# = 115 GeV, a correction AS" = 50 shifts the metastability bound 



by 1 GeV in m t . For this reason, our results are numerically close to those obtained in []T2 
(at zero temperature), where this correction was not included. Note that one could have 
expected a larger effect, because of the large value of the top Yukawa coupling g t at the 
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weak scale. However, the relevant quantity here is gt(l/R), with 1/R ~ 10 16 GeV, which 
is substantially smaller than g t (mt). 

The uncertainties in the determination of AS" turn out to be negligible with respect 
to the uncertainties involved in the determination of A(/x) at high scales. We estimate the 
latter to generate an error 6m t ~ 2 GeV at fixed m# and a s , both in the stability and 
metastability bounds, as shown in Eqs. ( |5.2j ) and (|5.3|). 

Some comments are in order: 

• Vacuum decay can also be catalysed by collisions of cosmic rays. However, it has 
been shown [[23[ that the tunnelling rate induced by such processes is negligible with 



respect to the probability of quantum tunnelling. 

Thermal tunnelling gives stronger bounds than quantum tunnelling only under the 
assumption that the temperature of the universe has been above 10 s — 10 9 GeV (for 
m H = 115 GeV). 

We have assumed the validity of the Standard Model up to the Planck scale. Ex- 
periments at Tevatron run IIB can reduce the error on m t down to ±2 GeV [E5 and 



possibly discover the Higgs and measure its mass if m# < 130 GeV. Some extension 
of the SM will become necessary, should m# and m t be found in the excluded region. 
All concrete modifications invoked to cure this problem give a computable correc- 
tion to the squared Higgs mass; for example, a scalar with mass fh and a coupling A 
to the Higgs gives 5m 2 H ~ \fh 2 / (An) 2 . Since A cannot be too small, by naturalness 
arguments we expect the scale fh of new physics to be in the electroweak range. 

All observed neutrino anomalies could be explained in terms of neutrino oscillations 



without affecting our result f24|. 



6 Conclusions 

All recent LEP1, LEP2, SLD and Tevatron data are compatible with the Standard Model 
with m t ~ 175 GeV and a light Higgs, maybe m# « 115 GeV 0. It is well known that 
the Standard Model is affected by the so-called hierarchy problem, which manifests itself 
through uncomputable quadratically divergent corrections to the Higgs squared mass. 
The extensions of the theory proposed to cure this potential problem are likely to involve 
the presence of new phenomena not far above the electroweak scale. However, a definite 
solution of the naturalness problem is not yet established, and it is interesting to consider 
the possibility that it be solved by some unknown mechanism that takes place around the 
Planck scale. Even within this assumption, a possible phenomenological problem affects 
the SM, namely the possibility that the scalar potential become unbounded from below at 
large field values, below the Planck scale. For m# = 115 GeV such instability is present if 
the top mass is larger than (166 ± 2) GeV, i.e. only about two standard deviations below 
the central value of the Tevatron data. 
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The instability of the SM vacuum does not contradict any experimental observation, 
provided its lifetime is longer than the age of the universe. In semi-classical approximation, 
this happens if the running quartic Higgs coupling A is larger than about —0.05. In this 
paper we have presented a more accurate assessment of this bound by performing a full 
next-to-leading order computation of the tunnelling rate of the metastable vacuum. In 
particular, we have computed the one-loop corrections to the bounce action. Our result is 
summarized in Eq. ( |3.3|) and plotted in Fig. f|. Next-to-leading corrections turn out to be 
numerically small, with an appropriate choice of the renormalization scale (/i = 1/R), but 
they allow fixing the various ambiguities of the semi-classical calculation, thus reducing 
the overall uncertainty of the result. We find that, for m# = 115 GeV (that is, just above 
the present exclusion limit), the instability is dangerous only for m t > (175±2) GeV; this 
result does not depend on the unknown new phyiscs at the Planck scale. A determination 
of m t at this level of accuracy would therefore be extremely interesting in this respect, 
and will probably be achieved at Tevatron run II. 
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A Spherical harmonics 

The spherical harmonics Yj mm i{6) are eigenfunctions of L 2 = L^L^, where 

i 

L^ u = (Xfidv - Xud/j) , (A.l) 

with eigenvalues 4j'(j+l); the indices m,m' range between — j'andj (withj = 0,1/2,1,...) 
giving the multiplicity (2j + l) 2 . This multiplicity can be understood by noting that the 
group of rotations in four dimensions is 0(4) ~ SU(2) (g) SU(2), and its scalar represen- 



tations correspond to SU{2) ® SU(2) representations (/, k) with / = k p6fl . The indices 
m, ml will be omitted in the following. 

Spin-1/2 representations are obtained for I — k ± 1/2. The index J that appears in 
Eq. ( |4.34| ) is related to / and k through 

J = 2k + - for / = k + -, k = 0,1/2,1,... ; (A.2) 

2 2 

J=-2k-- for / = k--, k = 1/2, 1, . . . . (A.3) 

2 2 

Therefore, the multiplicity of the representation labelled by J is (2/c+l)(2/ + l) = J 2 — 1/4. 
The action of the operator ^)on the (/, k) spinor representations can be found in Ref. [IHJ. 
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In the case of vector representations, it is straightforward to work out the action of d 2 
on A^x) as in Eq. Q4.44Q . Using Eq. (fO|) we find 



d 2 ai j x fJ Y j (9) 
Therefore, we have 



a" + ^ a 12 + 3 
a' 2j L 2 -l 

J 



xMffl + au-dpYtf) (A.4) 
2L 2 

rdpYjW + av—xMe) (A.5) 



(X;,: ~\~ 3 Clin n 



r 



" % 3 r 2 



-lip pa 



tfxAp&Ap = I r 3 dr <Hj(r)D$a kj 



(A.6) 



(A.7) 



where, with an appropriate normalization of the two vectors 



d 2 3d 4j(j + l) 



D {j) = + 

dr 2 r dr 







r 2 






£ 





(Ai 



B Fourier transform of the bounce 

The four-dimensional Fourier transform of a function /, depending only on the radial 
coordinate r, is given by 



f(p 2 ) = / d A xe ip ^f(r) = / drr 2 f(r).h(pr) 



p Jo 



where p 2 = p^p^ and 



In the case of the bounce we find 

Hp 2 ) r 
h\ v 2 ) 



JJpr) = ?L r d9 sm 2 6e iprcose . 

TV JO 

8V2tt 2 R 2 K 1 { P R) 



p Jo 

47T 2 roo 



drr h(r)Ji(pr) 



pj | A 



drr h (r)Ji(pr) 



l6n 2 R 2 



■KoipR) 



p Jo |A 
The functions K n (x) and J n (x) are Bessel functions defined as in Mathematica [27] 



(B.l) 
(B.2) 

(B.3) 
(B.4) 
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